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Abstract 

The magnetic order resulting from an indirect exchange between magnetic moments provided by 
spin-polarized hole gas in the metallic phase of a Gai_j;Mna;As double layer structure is studied via 
Monte Carlo simulation. The coupling mechanism involves a perturbative calculation in second 
order of the interaction between the magnetic moments and carriers (holes). We take into account 
a possible polarization of the hole gas due to the existence of an average magnetization in the 
magnetic layers, establishing, in this way, a self-consistency between the magnetic order and the 
electronic structure. That interaction leads to an internal ferromagnetic order inside each layer, 
and a parallel arrangement between their magnetizations, even in the case of thin layers. This fact 
is analyzed in terms of the inter- and intra-layer interactions. 
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I. INTRODUCTION 



During the last decade a new interest arose in the study of the magnetic order in layered 
materials. This area includes the study of magnetic semiconductor pseudo-binary alloys like 
Ai.^Mj-B, where M stands for a magnetic ion. These alloys are called Diluted Magnetic 
Semiconductors (DMS).0i A few years ago some groupg^l&DQ'QlaB succeeded in producing 
homogeneous samples of Gai_i.Mn^.As alloys with x up to 7% using low temperature (200 — 
300° C) Molecular Beam Epitaxy (MBE) techniques. Mn is a transition metal having its 3d 
level half filled with five electrons, in such a way that it carries a spin of 5h/2, according to 
the Hund's rule. In Gai_a;Mnj.As alloys a substitutional Mn acts as an acceptor (it binds one 
hole), and at the same time it carries a localized magnetic moment, due to its 3d shell. For 
X = 0.053, the alloy is a metallic ferromagnet,! the Curie- Weiss temperature is llOK, and 
the free hole concentration is near 1 — 2 x lO^'' cm~^. The ferromagnetic order in the metallic 
phase is understood as resulting from the indirect exchange between the Mn ions due to the 
local spin polarization in the hole gas.0 At present a vast literature exists about different 
possible mechanisms for that indirect exchange The free hole concentration in the 
metallic phase is a fraction (10-20%) of the total concentration of Mn. This is understood 
as due to the presence of As anti-sites and interstitial Mn. The magnetic ordering resulting 
from indirect exchange via spin-polarized free carriers implies in a spin coherence length 
larger than the average distance between localized magnetic moments. The possibility of 
having a DMS based on GaAs opens a wide range of potential applications as, for instance, 
in integrated magneto-optoelectronic devices. 

In this work we extended a confinement-adapted Ruderman-Kittel-Kasuya-Yosida 
(RKKY)0 mechanism to study the magnetic order resulting from the indirect exchange 
between magnetic moments in a GaAs/Gai_2,.Mn2:As nanostructure with two DMS layers. 
A temperature dependent Monte Carlo (MC) simulation is performed to determine the re- 
sulting magnetic phases. This article is organized as follows. In Sec. |I| we present the 
calculation of the indirect exchange for a confined spin-polarized Fermi gas in a semiconduc- 
tor heterostructure. In Sec. |IT| the self-consistent calculation of the heavy-hole single band 



electronic structure is presented in some detail. In Sec. |^ the Monte Carlo simulation is 
performed to determine the resulting magnetic phases and the relevant properties. In Sec. [V| 
we present our conclusions. Our calculations reveal that a ferromagnetic phase exists even 
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in thin layers in the structure containing two DMS layers. Considering the intra-layer and 
the inter-layer interactions as two independent mechanisms, we conclude that the intra-layer 
interaction is dominant in what concerns the internal magnetic order of each layer, while 
the inter-layer interaction alone, which determines the relative orientation of the two aver- 
age magnetizations, is responsible for a transition temperature four times lower than that 
resulting from the two mechanism acting together. The roles of the carrier concentration 
and of the cutoff parameter are also analyzed. 

II. INDIRECT EXCHANGE VIA SPIN-POLARIZED FERMI GAS 

The indirect exchange between localized magnetic moments in quasi-two-dimensional 
structures mediated by a Fermi gas has been addressed several times.0'0'0Eilli0illii Basi- 
cally, it deals with a confined electron (or hole) gas being locally spin-polarized by magnetic 
moments distributed in a layer. 

The interaction potential between the Fermi gas and the set of localized magnetic mo- 
ments is well described by the Kondo-like exchange term: 

i^ex = -IY.S^■ S(f)5(f - R,), (1) 

i 

where the localized spin of the Mn ion Si at position Ri will be treated as a classical variable, 
and s{f) is the spin operator of the carrier at position f; I is the sp — d interaction. 
If '?/'o-(r) and ipl.{r) describe the fermion field operator for spin a, then 

= l{^{r)Mr) - ^l(r1^i(f)), (2) 
s+(f) = V^|(f)V^l(f), (3) 

s-(f) = ^{(f)^/'|(f), (4) 

with the usual definitions of s"*" = -|- isy, and s~ = — isy. In terms of field operators, 
the (Kondo) exchange term in the Hamiltonian is 

^ex = — T.{St[i^W^)i^^i^i) - ^I(n)V^i(r1)] + St^imm + S.ij\{n)ijm}. (5) 
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We consider holes in a semiconductor heterostructure that are confined in the growth 
direction, assumed to be the z-axis. The total Hamiltonian, H = Ho + Hex includes in Hq the 
kinetic part, the confinement potential and the Hartree as well as exchange and correlation 
terms.ii Neglecting scattering from impurities, holes in the effective mass approximation are 
free particles in the {x,y) plane, i.e., in the plane parallel to the layer interfaces. Their field 
operator can be written as 

Mr) = ^ E e'%nA^)v.c^rk,^, (6) 

where A is the normalization area. A; is a wave vector in the {x, y) plane, r/o- is the spin 
tensor for the polarization a, (f)n,a{z) is the envelope function which describes the motion of 
the fermion in the z-direction, and c^j:^ is the fermion annihilation operator for the state 
(n, fc, cr). p represents a vector in the {x,y) plane. 

The confined RKKY indirect exchange is a second-order perturbative treatment. It gives 
two terms for the correction to the ground state energy of the system formed by the set 
of (classical) localized moments and the Fermi gas.0 The first one is the self-energy term, 
which corrects the site energy. The second is the interaction between the localized magnetic 
moments, the RKKY exchange, allowing the interaction to be written as a Heisenberg 
Hamiltonian: 

SE'^'^ = -J2J^A■s,. (7) 

The RKKY approach can be extended to spin-polarized states. This is, in one sense, going 
beyond the second-order perturbation. In that case no simple spin-spin scalar product is ob- 
tained, as in Eq. (|^), since the spin polarization breaks the rotational symmetry establishing 
a preferential direction, namely that of the average magnetization. 

In the second-order perturbation, the correction to the energy which is bilinear in the 
magnetic ion dipole moment is obtained from 

=- 1^ p — zr~^o • v8j 

t',F' ^F',t' ^F,t 

Here, the state \F,t > is understood as the direct product of the state of the system of 
classical localized moments and that of the quantum spin-polarized gas. Substituting Eq. @ 
into Eq. (R) we obtain, after a lengthy calculation, the effective Hamiltonian which includes 
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explicitly the spin-flip terms: 

Heff = - E(CJ/ + CiDS^S^ + {CJI + Ci]){S^SJ + SfS]), (9) 



where the spin-flip term is given by: 



Ci/ = - E E E(:^)' exp[zq(R, - R,)] x 



bnizj)e{Ep - e„_k)^(e„, k+q - Ep). (10) 

Notice that the arrows refer to spins up or down for holes, not the orientation of the localized 
magnetic moments at the Mn sites. Similar definitions hold for the other C"'"'' . Eq. (p^O]) can 
be expressed in terms of the Lindhard function:i3il 



^ ^n',k+q ^n,k 



(11) 



and its real-space Fourier transform: 

x"'^'(R.,) = EexpHg ■ R.,]x"'"'(g). (12) 

Substituting Eq. (|1T]) in Eq. (^) we obtain 

= - E E (^)'C'(^.)0n(;^.)0:(%-)0n'(%)x"'"'(R.,), (13) 

where the summations on n and nf are restricted to those sub-bands with the proper spin 
polarization. 

The intra-subband real space Lindhard function x"'"(Rjj) has been derived by several 
authors:0Sllilii 

X"'"(^i.) = ~^^kP'[MkPR,,)No{k^;^^Ri,) + MkPR,,)m{kPR,^)l (14) 

where is the transversal (parallel to the interfaces) effective mass, and k^p^ refers to the 
Fermi wavevector of subband n. 

The inter-subband terms cannot be expressed in a closed form.0 They must be obtained 
numerically from the integral 

X"'" (Ri,) = - / dqqF^,^,{q)MqRi,), (15) 
Jo 

where 



i^nx(g) = ^(l-%^)[l-\l-( m^' + ^n',n-2qkPme.^,-E,), (16) 

with Ar,',n = 2m* ■ {En' - E^)/h\ 



III. SPIN-POLARIZED ELECTRONIC STRUCTURE FOR MAGNETIC MULTI- 
LAYERS 



In order to obtain the spin-polarized electronic structure for holes given an average mag- 
netization of the Mn ions, we solve self-consistently the heavy hole single band Schrodinger 
equation in the reciprocal space. The hole system is supposed to be homogeneous in the 
xy plane, so the Hartree term depends only on the coordinate z, Unir) = Uh{z). For the 
purpose of obtaining the electronic structure we treat the magnetic interaction as being 
due to an uniform magnetization in the DMS layers. If a net magnetization exists, it will 
polarize the hole gas. This problem is solved self-consistently by a secular matrix equation 
in the reciprocal space. The method would be exact were not for cutting the matrix size. 
The advantage is that it provides spin-polarized eigenvalues and eigenf unctions with high 
accuracy, not only for bound states, but also for a high number of scattering states. For 
each spin, we define the wavefunction Fourier Transform (FT): 

i^air) = J £qexp{iq.r)'ip^{q). (17) 
The hole eigenstates will be obtained by discretizing the integrals on q appearing in 

J d^r^;{r){H - E)%l)„{f) = 0. (18) 

When integrating the magnetic term in the Hamiltonian over r, we assumed the magnetic 
impurities to be uniformly distributed in each one of the Gai_j.Mn^As DMS layers, all of 
them having the same thermal average magnetization. This treatment includes not only 
the ferromagnetic phase but also phases where a partial magnetization is observed, as the 
"canted-spin" phases.ll2B Therefore, taking an homogenous concentration of Ni magnetic 
impurities in the j-th layer, we have: 

-I rfVexp[i(g-g').rl^^s(f).S^(^i)5(f-i?i) ~ 

Ni 

-I-J2<M >, Eexp[^(g - (^).R^] = 

-/|^ < M ^ J d'R,exp[t{q- (f).R,] = 
j 
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where a = ±1 for spin parallel (upper sign) or anti-parallel (lower sign) to the magnetization, 
and y- = a;A^o is the impurity density, in terms of Nq, the density of the cation ions. The 
sp-d exchange constant for holes is usually written as (3. Here we used Nq/S = — 1.2eV.izl 
^DMS is the integral performed on the z -coordinate inside the j-th DMS layer: 

FDMsil) = ^ I dze^^[iq.z\. (20) 
In Jj 



Comparing Eqs. ([19| ) and (l20|), the latter being just the Fourier transform of an unit 
barrier function within the range of the j-th DMS layer, we see that the thermal average 
magnetization < M >j polarizes the hole gas by introducing additional effective confining 
potentials given by 

= -No/3x^ E < ^ 9A^)^ (21) 
^ j 

where gj{z) = 1 if 2; lies inside the j-th layer, and gj{z) = otherwise. Next, we define 

Ueffiq) = ^J dzexp[tq.z]mz) + V;JJ^iz) + Uniz)], (22) 

where Uc{z) and Uh{z) are the confining and the Hartree potentials, respectively. The 
eigenvalues and eigenfunctions at the bottom of the 2-D subbands may be obtained by 
solving the secular equation for each spin polarization: 



det 



2m* 



E 



g;) + t^e//(g.-g;)| = 0. 



5(g,-g;) + [/e//(g.-g;) =0. (23) 



IV. MONTE CARLO SIMULATION: INTERLAYER INTERACTION AND MAG- 
NETIC ORDERING 

In the present work we have focused our attention on metallic systems with the Mn 
concentration x = 5%, and we have neglected possible superexchange (anti-)ferromagnetic 
interaction between the nearest neighbors and the next-nearest neighbors pairs. We have 
performed extensive Monte Carlo (MC) simulations in order to determine the possible mag- 
netic order occurring in a system containing a double layer of Gai_j;Mna;As. Classical spins 
Si, representing the localized magnetic moments of the Mn ions, are randomly distributed 
on the cation sites with concentration x. They are assumed to interact through the exchange 
Hamiltonian defined by Eq. (^. 

The interaction derived in Sec. |I| is assumed to be effective within a cutoff radius which 
corresponds to the carrier spin coherence length. This cutoff is taken tentatively as Rc = 



2.5a, and Rc = 2a, where a is the fee lattice parameter of GaAs. These values correspond to 
five and four monolayers (ML), respectively. This choice makes the smallest value assumed 
for Rc nearly equal to the hole mean free path estimated from bulk transport measurements. I 
The electronic structure was calculated for a system of two DMS layers grown inside a 
GaAs quantum well. We have studied two different DMS layer widths, 10 and 20 A. The 
GaAs spacer separating the DMS layers has been considered as having widths varying from 
5 A to 60 A. The whole system consists of a 200A quantum well with infinite barriers in the 
z-direction, in the center of which the two GaMnAs layers and a GaAs spacer are disposed, 
as shown in Fig. ^ The eigenstates, obtained according to Sec. were calculated for 
net magnetizations = 0.0, 0.1, 0.2, . . . 1.0, for each one of the geometries used. Notice 
that < M > determines the effective magnetic potential in Eq. and, in consequence, 
the whole electronic structure. A typical result for the electric charge density, and the spin 
polarization densities is shown in Fig. Ic. We see that holes are attracted to the region of the 
DMS layers by the negatively charged Mn ions. Because Nf)(3 is negative, however, the sp — d 
interaction (taken into account as the effective magnetic potential) favors the occurrence of 
anti-parallel spin subbands at lower energies. Several subbands are occupied, due to the 
high concentration of free carriers, most of them corresponding to the anti-parallel spin 
polarization. Therefore, the total interaction, i.e., the confinement, the Coulomb interaction 
(Hartree and correlation terms) plus the sp — d interaction results in an inhomogeneous 
distribution of spin polarization and spin charge density. Holes with parallel spins, at lower 
density than that of parallel spins, concentrate on the DMS layers, while the latter, at 
higher density, spreads in a much wider region. The inhomogeneity of the spin polarization 
density influences the indirect exchange between localized magnetic moments, and this fact 
is taken into account explicitly in the present calculation. For each value of < M > the 
exchange interaction terms in Eq. (^ are calculated and tabulated for all possible magnetic 
ions distances Rij in the (x, y) plane according to the zincblende structure within the cutoff 
radius. 

The MC calculation is performed in a finite box, whose axes are parallel to the [100] 
directions. Its dimensions are = Ly, and = 2L, where L is the width of each DMS 
layer. Periodic boundary conditions are imposed in the (x, y) plane. The lateral dimensions 
are adjusted in such a way that the total number Ng of Mn sites is about 4000, for all samples 
with different L^. The initial spin orientations of the Mn ions are randomly assigned in the 
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two DMS layers which are separated by a distance the width of the GaAs spacer. At a 
given temperature, the energy of the system due to exchange interaction is calculated, and 
the equilibrium state for a given temperature is sought by changing the individual vector 
spin orientation according to the Metropolis algorithm.!! A slow cooling stepwise process 
is accomplished, starting from above the transition temperature T^, and making sure that 
the thermal equilibrium is reached at every temperature. Every time the net magnetization 
increases to reach a value of n x 0.1, n an integer varying form to 10, the tabulated 
values of the exchange interaction terms are substituted by those corresponding to the new 
reference < M >. The tabulated values of the exchange interaction terms were calculated 
from the hole states (eigenfunctions, eigenvalues and Fermi wavevectors) resulting from self- 
consistent calculation described in Sec. |IT| for this specific thermal average magnetization 
< M >. Then, the resulting spin configuration for the Mn ions is taken as the starting 
configuration for the next step at a lower temperature. The Monte Carlo procedure adopted 
in this way takes account of the changes on the hole gas polarization due to the presence 
of the established order among the localized magnetic moments in the Mn ions. Thus 
the magnetization is determined by a fully self-consistent process, but the effects of the 
fluctuations of the magnetization on the electronic structure are neglected in this process. 

For each temperature, the thermal average magnetization < M >, and the Edwards- 
Anderson (EA) order parameter q are calculated.!! The latter is defined as 



where a = x, y and z. In order to avoid spurious results on the thermal average values of q 
over a large time interval t, a summation on t' is performed starting from a time to, when 
the system already reached the thermal equilibrium. This is true for other quantities such 
as energy, magnetic susceptibility, etc. 

The MC simulations have been performed for several samples as indicated on Table I. 
The ratio between the carrier concentration and the Mn concentration is the same for the 
two DMS layers, and is chosen to reproduce hole concentrations which in bulk would be 
p = 1 X 10^° cm~^, and p = 2 x 10'^^ cm^^. In Fig. the normalized magnetization as a 
function of temperature is shown for samples ^^1 to ^6, with the DMS layer width of 10 
A, bulk hole concentration p = 1 x 10^° cm~'^, and Rc = 5 monolayers (ML). All these 




(24) 
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samples show a paramagnetic (P) to a ferromagnetic (FM) phase transition. The transition 
temperatures decrease with the increase of the spacer thickness, pointing to the importance 
of the inter-layer interaction. In Fig. |^ all parameters are kept the same as in the previous 
figure, but p = 2x 10^° cm~'^. The increase on the hole concentration has two main influences 
on the magnetization curves: (i) an increase of and (ii) the possibility of the occurrence 
of a canted spin phase (C), with a partial ferromagnetic order .0 This canted spin phase 
can be explained on the basis of the oscillatory behavior of the exchange interaction. The 
increase on the hole concentration produces an increase of the Fermi wave number, and 
anti-ferromagnetic (AF) interactions are likely to occur, according to the oscillatory nature 
of the (confined) indirect interaction. The Curie temperatures for samples 7^ 01 to 7^ 12 
are obtained in the range from 20 K to 49 K. 

In Fig. ^ the results are shown for two 20 A DMS layers. The hole concentration is 
j9 = 2 X 10^° cm~^, and the cutoff length is 5 ML. decreases with the interlayer separation, 
as before. However, compared to the 10 A system, a higher is observed in the present 
case. In Fig. |^ the cutoff radius used was 4 ML. We observe here the same behavior as 
before, but with higher T^, and also with higher values of < M >. The presence of effective 
AF interactions in these samples, due to the oscillatory behavior of the exchange interaction, 
may produce canted spin phases, rather than FM phases. This results from competing ferro- 
and antiferromagnetic interactions. It is the case of the so-called frustration. Although these 
are ordered phases, as can be observed in the EA parameters on Fig. ^, the Mn magnetic 
moments are not all parallel and the magnetizations at zero temperature are in most samples 
partial, lying on the range of 60% to 100%. The Curie temperatures for samples 7^ 13 to 
7^ 24 are obtained in the range from 35 to 64 K. 

In order to better understand the importance of the spin-polarized indirect exchange on 
the magnetic order and on the transition temperature of these layered systems, we have 
compared our results with a simulation performed with an unpolarized exchange. This is 
done by using the tabulated values of the exchange for < M >= in the whole cooling 
stepwise MC simulation. For this comparison we have chosen a sample with the same 
characteristics of sample 7^ 1, i.e., two 10 A DMS layer and a GaAs spacer of 5 A. The result 
is shown in Fig. |^. The presence of polarization more than doubles the Curie temperature. 
This occurs because the effective magnetic potential resulting from the finite thermal average 
magnetization, and the consequent polarization of the hole gas favor a stronger concentration 



10 



of charge and spin on the DMS layers, as can be seen in Fig. |l|, providing a stronger exchange 
than that for the unpolarized system. 

A calculation was performed to estimate how the interlayer interaction is affected by the 
the width of the spacer, and how this effect appears in the Curie temperature. To this end, 
the exchange interaction inside each DMS layer was "turned off" by making it artificially 
equal to zero. Thus, only the Mn ions belonging to two different DMS layers are interactive. 
The MC simulations are then performed on the same basis as described before. The results 
are presented in Table II, for the two DMS layers with 10 A width separated by 5, 10, 
15, 20, 40, and 60 A, and for two different hole concentrations, 1 x 10^'^ cm~^ in Fig. |^ 
(samples #19 to #24), and 2 x 10^° cm^^ in Fig. | (samples #25 to #30). These results 
show that there is an effective interlayer coupling able to produce a FM phase. The fact 
that Tc coming out exclusively from that interaction is about four times lower than that 
for the complete interaction reveals that the main contribution to the magnetic ordering 
within each layer is due to de intralayer coupling. The interlayer interaction acts mainly 
to establish the ferromagnetic arrangement between the magnetization of the two layers. 
Notice, however, that the resulting magnetic order affects and is affected by the distribution 
of charge and spin of the free carriers in the structure. We have performed an exponential 
fitting of the transition temperature as a function of the spacer thickness for samples with 
L = lOA , taking into account the complete interaction and also solely with the interlayer 
interaction. This is shown in Fig. TU. Notice that the exponential fitting works well in 
both cases, putting into evidence the importance of the intra-layer interaction. A similar 
calculation was performed to explore the dependence of the Curie temperature on the DMS 
layer width, for a fixed spacer thickness {d = 10 A ). The results appear in Fig. |ll], for 
L = 10,20,40 and 50A. The Curie temperature increases with the DMS width approaching 
the fully filled T„ saturation exponential function. 

Up to now the simulations in this work were performed in a single Monte Carlo sample (for 
each geometry chosen, given by the set of parameters L, d and p) with a specific distribution 
of the sites occupied by the magnetic ions. Next, the calculation is performed taking into 
account a configurational average on the initial distribution of the orientation of individual 
magnetic moments Si on the Mn ions, as well as on the sites Ri where these ions are 
localized. The positions of the Mn ions identify a MC sample. The parameters are chosen 
as those of sample #01 in Table I. The site configurational average is performed, in each 
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case, over three initial configurations, namely a,b and c (the MC samples). Although the 
results do not change significantly, as it can be seen from Fig. |T2| which shows the average 
magnetization versus temperature, it is important to observe general behaviors which appear 
in the process. For instance. Fig. O shows the susceptibility versus temperature, and we can 
see that, apart the critical fiuctuations to be expected at the transition temperature, which 
is perfectly defined in that figure, the fiuctuations are more important at low temperatures 



than at high temperatures above Tc- Another interesting behavior appears in Fig. 
showing the magnetic specific heat (Cy) versus temperature. Again, the fiuctuations are 
more important at low temperatures. Besides, the maximum of Cy occurs at a temperature 
which is lower than that of the maximum of the susceptibility, what is reminiscent of the 
spin-glass behavior, and can be understood as a consequence of the competing FM and AF 
interactions characteristics of the indirect exchange. 



V. DISCUSSIONS AND FINAL COMMENTS 



This Monte Carlo simulations shows that this system with two thin GaMnAs layers inside 
a GaAs matrix presents a ferromagnetic, or at least a canted magnetic phase ordering. The 
input parameters, namely the hole concentration, the DMS layer width, the DMS layer 
separation and the interaction cutoff, all influence the magnetic phase of these systems. The 
increase of hole concentration favors the canted phase due to the intrinsic frustration of the 
indirect exchange interaction. In the reverse direction, the elimination of the AF interactions 
produces FM samples and an increase of the Tc as already observed in a previous work.0 The 
increase of the DMS layer width will increase the number of interacting magnetic moments 
raising also T^. As expected, T^, will be lowered with the increase on the spacer thickness, 
since as the separation increases the interlayer coupling become weaker and weaker until it 
reaches a limiting value. This limit is flnite due to the conflnement by the inflnite barriers 
at the boundaries of the structure. We have quantifled the influence of this parameter on 
the magnetization of the system. Despite the fact that this interlayer interaction is four 
times weaker than that of the complete system, it is important to notice that it plays an 
important role in the existence of a FM phase in multilayers, and in digital superlattices. 
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FIG. 1: (a) Structure model used in the present calculation: two DMS layers inside a GaAs 
quantum, (b) The effective magnetic potential (schematic), (c) Density distribution: charge(solid 
line), spin polarization(dashed), anti-parallel spin (dashed-dotted), and parallel spin (dotted). 

FIG. 2: Normalized magnetization vs temperature for samples ^^01 to ^^06 indicated in Table |. 

FIG. 4: Normalized magnetization vs temperature for samples #13 to #18 indicated in Table |. 

FIG. 5: Normalized magnetization vs temperature for samples #19 to #24 indicated in Table |. 

FIG. 6: Edward-Anderson order parameter vs temperature for samples #13 to #18 indicated in 
Table |. 

FIG. 7: Comparison for the polarized (solid line) and unpolarized (doted line) system. The char- 
acteristics of the sample are the same as sample #01, indicated in Table |[ 

FIG. 8: Normalized magnetization vs temperature for samples #25 to #30 indicated in Table || 
taking into account only the interlayer interaction. 

FIG. 9: Normalized magnetization vs temperature for samples #31 to #36 indicated in Table || 
taking into account only the interlayer interaction. 

FIG. 10: Curie temperature vs spacer thickness for L = lOA , and p = 1. x 10^^ cm~^ for the 
complete interaction (samples #01 to #06 indicated in Table ||) and for the interlayer interaction 
(samples #25 to #30 indicated in Table ||). An exponential fitting y{x) = yo + ^exp(— j;/A) 
is shown, with the resulting parameters: yo = 20. IK, A = 31. 2K, A = 9.8A for the complete 
interaction; yo = 5.9K, A = 15. 4K, A = 6.5A for the interlayer interaction. 

FIG. 11: Curie temperature vs DMS layer thickness for d = lOA , and p = 2. x 10^^ cm^^ An 
exponential fitting y{x) = yo{l — exp(— x/A)) is shown, with the resulting parameters: yo = 85. 4K, 
and A = 15.21 . 

FIG. 3: Same as Fig 2 for samples #07 to #12 indicated in Table |. 
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TABLE I: Sample characteristics: L is the width of each DMS layer; d is the GaAs spacer width; 
the ratio of the carrier concentration to the Mn concentration, written in terms of a bulk hole 
concentration, p ; Tc is the transition temperature for phases: FM: ferromagnetic, C: canted spin; 
Rc is the cutoff radius of the exchange interaction in number of monolayers. 



sample L (A) d (A) r Rc phase Tc (K) 



#01 


10 


05 


1 


5 


FM 


38 


#02 


10 


10 


1 


5 


FM 


31 


#03 


10 


15 


1 


5 


FM 


28 


#04 


10 


20 


1 


5 


FM 


23 


#05 


10 


40 


1 


5 


FM 


21 


#06 


10 


60 


1 


5 


FM 


20 


#07 


10 


05 


2 


5 


c 


49 


#08 


10 


10 


2 


5 


FM 


39 


#09 


10 


15 


2 


5 


FM 


35 


#10 


10 


20 


2 


5 


FM 


34 


#11 


10 


40 


2 


5 


FM 


32 


#12 


10 


60 


2 


5 


FM 


22 


#13 


20 


05 


2 


5 


c 


62 


#14 


20 


10 


2 


5 


C 


63 


#15 


20 


15 


2 


5 


c 


59 


#15 


on 


on 


o 
z 


5 


U 


Ob) 


#17 


20 


40 


2 


5 


c 


53 


#18 


20 


60 


2 


5 


c 


51 


#19 


20 


05 


2 


4 


FM 


64 


#20 


20 


10 


2 


4 


FM 


56 


#21 


20 


15 


2 


4 


FM 


52 


#22 


20 


20 


2 


4 


FM 


49 


#23 


10 


40 


2 


4 


FM 


35 


#24 


10 


60 


2 


4 


FM 


45 
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TABLE II: Same as table 1, but considering only interlayer interaction. 

sample L (A) d (A) p (lO^o cm'^) phase (K) 



#25 


10 


05 




5 


FM 


13 


#26 

II 


10 


10 




5 


FM 


09 


#27 


10 


15 




5 


FM 


08 


#28 


10 


20 




5 


FM 


06 


#29 

II — <^ 


10 


40 




5 


FM 


06 


#30 

It 


10 


60 




5 


FM 


06 


#61 


1 n 


Uo 


o 
z 





r iVi 


1 K 

io 


#32 


10 


10 


2 


5 


FM 


08 


#33 


10 


15 


2 


5 


FM 


08 


#34 


10 


20 


2 


5 


FM 


08 


#35 


10 


40 


2 


5 


FM 


08 


#36 


10 


60 


2 


5 


FM 


07 



FIG. 12: Thermal average magnetization as a function of temperature for three different mag- 
netic moments configurations (open symbols) , and the corresponding configurational average (solid 
square), corresponding to the parameters of sample #01. 

FIG. 13: Magnetic susceptibility vs temperature, calculated from thermal fluctuations, for three 

different magnetic moments configurations (open symbols), and the corresponding configurational 
average (solid square), corresponding to the parameters of sample #01. 

FIG. 14: Magnetic specific heat vs temperature, calculated from thermal fluctuations, for three 
different magnetic moments configurations (open symbols) , and the corresponding configurational 
average (solid square), corresponding to the parameters of sample #01. 
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Fig. 01 Boselli et al. 



fig. 02 Boselli et al. 
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fig. 03 Boselli et al. 
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fig. 04 Boselli et al. 
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fig. 05 Boselli et al. 
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fig. 06 Boselli et al. 
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Fig. 07 Boselli, et al. 




fig. 08 Boselli et al. 
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fig. 09 Boselli et al. 
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